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Finite-temperature investigation of quarter filled ladder systems 

C. Gabriel a , E. Sherman a , T. C. Lang a , M. Aichhorn a , H. G. Evertz a >* 

a Institute for Theoretical and Computational Physics, Graz University of Technology, Austria 



c3 



c 
o 
o 



> 

O 

o 

^— » 
03 



o 
o 



X 



Abstract 

We investigate charge ordering in a quarter-filled ladder at finite temperature by determinantal Quantum Monte 
Carlo. The sign problem is moderate in a wide range of model parameters relevant for NaV^Os. The charge order 
parameter exhibits a crossover as a function of inverse temperature /3 on finite systems. Above a critical nearest 
neighbor Coulomb repulsion V c , the correlation length grows exponentially with /3, indicative of the ordered phase 
at /? = oo. We find a clear single-particle gap manifesting itself in a flat n(n) dependence at large nearest neighbor 
Coulomb repulsion V. 
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The inorganic ladder compound NaX^Os has at- 
tracted great attention in recent years. This interest 
was triggered by magnetic susceptibility measurements 
[1], which show a phase transition at T — 34 K into a 
low-temperature spin-gapped phase. This transition is 
accompanied by charge ordering, as observed in NMR 
measurements [2] , where the valence of the vanadium 
sites changes from V 4 ' 5 to V 4,5± ' 5 , with 8 the amount of 
charge disproportion. This transition has been studied 
theoretically by several techniques at T = [3] . 

On a microscopic level the system can be described 
by an extended Hubbard model 

H = ~ X] iy(ctc 3V +H.c.) + t/X>itn a 

(ij) • 

at quarter filling (n) = 0.5, with hopping matrix ele- 
ments tij = t x along the ladder and ttj — t y within a 
rung, and chemical potential fi. We state all energies 
in units of t y . These hopping parameters as well as the 
onsite Coulomb interaction can be extracted from first- 
principle calculations [4]. The hopping along chains 



t x — 0.5t y is weaker than along rungs. This strongly 
influences the physics of the ladder, for which a spin- 
gap seems to appear at tx^J-y [3]. 

We used t x = 0.5 and U = 8. Since the non-local 
Coulomb interaction V cannot be determined properly 
by band-structure calculations, we used V as a free pa- 
rameter of the Hamiltonian. The charge order param- 
eter is A r 2 n = ^ T^ e iQ(r *-^(^ - {n))( nj - (n» 



with Q = (n, 7r), which is unity for complete ordering. 

We performed grand canonical calculations by deter- 
minantal quantum Monte Carlo. These are often very 
difficult for doped systems because of a sign problem. 
Fortunately, the average sign is favorably large in the 
relevant parameter range of t x /t y — 0.5 and large V 
(Fig. 1). In the opposite case of isotropic t x — t y at 
small V, (sign) becomes very small. The charge or- 
der parameter exhibits similar behavior, but it is less 
strongly dependent on t x /t y . Charge order grows with 
increasing V. 

Fig. 2 shows the charge correlation length £ C c- At 
small interactions, V — 1.5 and 2.0, the correlation 
length seems to saturate, but for V — 2.5 and 3.0 it 
increases exponentially with j3, with a ^-dependent 
slope. This behavior is consistent with the ID Ising 
model in a transverse field (IMTF) [5] , which is equiv- 
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Fig. 1. Mean value of sign and order parameter A co as functions 
of t x /t y at V — 3, — 6, L — 16, (n) — 0.5, and as functions 
of V at = 6 (inset). 




Fig. 2. Logarithm of the charge correlation length £ cc as a 
function of for different interactions V (L — 32 for V — 1.5 
and 2; L = 44 for V = 2.5 and 3) 

alent to Eq. 1 in the limit of one spinless electron per 
rung. For large V , the transverse field goes to zero, and 
£imtf = | lntanh(/3)| _1 . This is exponential behavior 
with slope 2 at large /3, which the data in Fig. 2 appear 
to approach. There is long range order in the thermo- 
dynamic limit only at f3 — oo. For weaker interactions, 
V < 2t y , the correlation length £imtf remains finite 
even in the limit /? — » oo, showing a disordered phase at 
all temperatures. The results in Fig. 2 are nicely con- 
sistent with recent DMRG calculations [6] which show 
that at T = the system has a quantum phase transi- 
tion to an ordered phase at V c = 2.1(1). 

However, the behavior of finite size systems is differ- 
ent from the IMTF in the thermodynamic limit. As a 
function of inverse temperature, the charge order pa- 
rameter exhibits a crossover at large V (Fig. 3). For 
the single ladder, this crossover is a finite size effect. It 
appears since at some /3, the charge order correlation 
length will exceed the system size, resulting in appar- 
ent long range order. For a three-dimensional system 
of coupled ladders, this crossover can become a phase 
transition. The order parameter at V — 3 and different 




Fig. 3. Charge order parameter as a function of for V — 3.0 
and V — 2.5 with L — 32 rungs. 
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Fig. 4. Mean electron density as a function of the chemical 
potential [i at — 6, L— 16. 

P scales well as a function of £(L)/L. 

The onset of charge order at large V is most clearly 
visible in the single particle gap shown in Fig. 4. It 
manifests itself as a plateau in the n(/x) dependence 
where (n) = 0.5 remains constant in a region /j, m ln < 
A* < Mmax- At large V ', the upper boundary /^ max shifts 
with V as approximately 3V, which is the same value 
as in the atomic limit at full ordering. 
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